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Abstract 

We briefly review the hadronization pictures adopted in the LUND String Fragmenta- 
tion Model(LSFM), Webber Cluster Fragmentation Model(WCFM) and Quark Combination 
Model(QCM), respectively. Predictions of hadron multiplicity, baryon to meson ratios and baryon- 
antibaryon flavor correlations, especially related to heavy hadrons at Z° factory obtained by LSFM 
and QCM are reported. 
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I. INTRODUCTION 



The hadronization mechanism is an important but still unsolved problem up to now 
due to its nonperturbative nature. It is recognized that the hadronization mechanism is 
universal in all kinds of high energy reactions, e.g., e + e~ annihilation, and hadron(nuclear)- 
hadron (nuclear) collisions. Among these reactions, e + e~ annihilation at high energies, es- 
pecially at the Z° factory in the future, is best for studying the hadronization mechanism, 
since all the final hadrons come from primary ones, all of which are hadronization results. 
The e + e~ — > 7*/Z° — > h's process is generally divided into four phases (see Fig. [p. 

1. In the electro- weak phase, e + e~ pair converts into a primary quark pair qq via virtual 
photon or Z°. This phase is described by the electro- weak theory. 

2. The perturbative phase describes the radiation of gluons off the primary quarks, and 
the subsequent parton cascade due to gluon splitting into quarks and gluons, and the 
gluon radiation of secondary quarks. It is believed that perturbative QCD can describe 
this phase quantitatively. 

3. In the hadronization phase, the quarks and gluons interact among themselves and 
excite the vacuum in order to dress themselves into hadrons, that is, the confine- 
ment is 'realized'. Since this process belongs to the unsolved nonperturbative QCD, 
investigations employing various models will shed light on understanding this process. 

4. In the fourth phase, unstable hadrons decay. This phase is usually described by using 
experimental data. 

This paper focuses on step 3. The main method is to compare results of various hadroniza- 
tion models with the data, and the future Z° factory is very suitable for this purpose. The 
popular hadronization models at the market, Lund String Fragmentation Model(LSFM)[l|, 
Webber Cluster Fragmentation Model(WCFM) 2, 3], and Quark Combination Model(QCM), 
succeed in explaining a lot of experimental data in e + e~ — > h's and pp(p) — > h's processes 
by adjusting corresponding parameters. Recently, the baryon to meson ratioj^j, [s| and con- 
stituent quark number scaling of elliptic flow t^jfi] are measured at RHIC experiments, 



which do not favor the fragmentation model, while the QCM can explain these phenomena 

On n 

naturally [7H10J. QCM was first proposed by Annisovich and Bjorken et al. [llj. It was famous 
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e + e" -> y*/Z° -> qq parton hadronization particle decays 

shower 



FIG. 1: e + e annihilation to hadrons. 



for its simple picture and its successful prediction of the percentage of vector mesons. One 
of its great merits is that it treats the baryon and meson production in an uniform scheme, 
so it describes the baryon production naturally. However, for a long period QCM is regarded 
as being ruled out, because the prediction for baryon-antibaryon(I?I?) rapidity correlation 
in e + e~ annihilation by Cerny's Monte Carlo Program, which was alleged to be based on 
QCM, has great discrepancies with the experimental results of TASSO collaboration 12]. 
But early in 1987, the BB phase space correlation from the naive QCM scheme was ana- 
lyzed, which showed that there should not be such an inconsistency qualitatively. In the 
meantime, Quark Production Rule and Quark Combination Rule(the so-called 'Shandong 



Quark Combination Model 



SDQCM)') were developed. Then a serials of quantitative re- 



sults obtained by SDQCM JjJ- 15] confirmed that SDQCM can naturally explain the BB 
short range rapidity correlation together with Baryon to Meson ratio when the multi-parton 
fragmentation is included. In order to understand the hadronization phenomena especially 
in heavy ion collisions, it is necessary to study the hadronization mechanism in detail in 
e + e~ annihilation once Z° factory is available. On the other hand, the hadronization model 
serves as a bridge between the perturbative QCD and experiments, so it is a very important 
tool for studying, e.g., LHC physics. The properties of the light hadrons have been studied 
in our previous works. Here we focus on investigating the production of heavy hadrons (e.g., 
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Ac, At, B u , etc.) by LSFM and SDQCM. 

This paper is organized as follows: In section |TT1 we give a brief introduction to the 
popular hadronization models, i.e., LUND String Fragmentation Model, Webber Cluster 
Fragmentation Model and Quark Combination Model. Some numerical results are given in 
section III II Finally, a short summary and outlook close the study. 



II. HADRONIZATION MODEL 



A. LUND String Fragmentation Model 



String Fragmentation Model, first proposed by Artru and Mennesser in 1974[16|, has been 
developed by the theory physics group of Lund University since 1978, and corresponding 
Monte-Carlo programs (e.g., JETSET, PYTHIA, etc.) are written. By now PYTHIA is one 
of the most widely used generators describing the high-energy collisions. 

The hadronization of qq color-singlet is the simplest case for String Fragmentation. Lat- 
tice QCD supports a linear confinement potential between color charges, i.e., the energy 
stored in the color dipole field increases linearly with the separation between them. The 
assumption of linear confinement is the starting point for the String Model. As the q q move 
away in the opposite direction, the kinetic energy of the system changes into the potential 
energy of the color string (or color flux tube). When the potential increases to a certain 
extent the string will break by the production of a new quark pair, and the production 



possibility is given by the quantum mechanical tunnelling 



17H19] 



P q , m = e ~™yn = e -n m 2 / Ke -np 2 T /K ^ 

where my is the transverse mass of q'(q') quark, and k, the string constant, denotes the 
potential energy per unit length. Considering the assumption of no transverse excitation of 
the string, the px is locally compensated between q' and q' . From Eq. ([1]), one sees clearly 
that the probability of the different flavors is about w:d:s:c~l:l: 0.3 : 1CU 11 . So c, b 
quarks are not expected to be produced through fragmentation, but only in the perturbative 
QCD procedure. 

When the new quark pair q'q' are excited out in vacuum, the qq color string splits into 
qq' and q'q two color-singlets. If the invariant mass of either of the singlet system is large 
enough, further break will occur. The process ends when all string pieces become exactly 



hadrons on their mass-shells. This picture describes meson production naturally other than 
the production of baryons. More complex diquark and popcorn mechanisms have to be 
introduced to explain the baryon production. One prominent feature of LSFM is that 
energy momentum and flavour are conserved at each step of the fragmentation process. 

In e + e~ — > h's process at high energies, multi-parton states will be produced at the end 
of perturbative phase, so that the multi-parton fragmentation has to be taken into account. 
The corresponding procedures have been developed in ref. (2^]. The fragmentation picture 
of LUND model is shown in Fig. [2j 




String 1 String 2 



FIG. 2: string fragmentation 



B. Cluster Fragmentation Model 

The Cluster Fragmentation Model has been proposed by Wolfram^] in 1980. Webber 
Cluster Fragmentation Model(WCFM)[3] is the well-known example. It has three parts: 

• The formation of color-singlet cluster. Gluons split into quark pairs after parton 
shower: g — > qq. Adjacent quark and anti-quark from different gluons combine to 
form a color-singlet. 

• Cluster with large invariant mass fragments into smaller ones. 

• Little cluster decays into primary hadrons. i.e., Cluster — > hadronl + hadron2. 

Note that only two-body decay or fragmentation is adopted in the Webber cluster fragmen- 
tation model. The corresponding hadronization picture is local, universal and simple. Its 
fragmentation picture is shown in Fig. [3j 
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FIG. 3: cluster fragmentation 



C. Quark Combination Model 



In the framework of SDQCM, Quark Production Rule(QPR), Quark Combination Rule 
(QCR) are adopted to describe the hadronization in a color singlet system, and then a 
simple Longitudinal Phase Space Approximation (LPSA) is used to obtain the momentum 
distribution for primary hadrons in its own system. Finally this hadronization scheme is 



extended to the multi-parton states. E 
relevant equations(for detail, see ref. [15 



ere we briefly introduce the SDQCM and list the 



2l|). 



1. QPR and QCR for qq system 

In a color singlet system formed by qq, N pairs of quarks can be produced by vacuum 
excitation via strong interaction. We assume that A" satisfies Poisson Distribution: 

n<N>,N-l)= < *> N iy \ -< N > (2) 

where < A^ > is the average number of those quark pairs. According to QPR, < A^ > is 
given by 

< N >= ya 2 + (5{W -M q - Mg + 2m) - a - 1, a = (3rh - - (3) 

where W is the invariant mass of the system, (3 is a free parameter, fh is the average mass 
of newborn quarks, and M q and Mq are the masses of endpoint quark and anti-quark. Thus 
we have A^ pairs of quarks according to eqs. (121), ([3]) (containing one primary quark pair). 

When describing how quarks and antiquarks form hadrons, we find that all kinds of 
hadronization models satisfy the near correlation in rapidity more or less. Since there is 
no deep understanding of the significance and the role of this, the near rapidity correlation 



has not been used sufficiently In ref. 



13], we have shown that the nearest correlation 



in rapidity is in agreement with the fundamental requirements of QCD, and determines 
QCR completely. The rule guarantees that the combination of quarks across more than two 
rapidity gaps never emerges and that N quarks and N antiquarks are exactly exhausted, 



which guarantees the unitarity (see ref. 22]). Given that the quarks and antiquarks are 



stochastically arranged in rapidity space, each order can occur with the same probability. 

Then the probability distribution for N quarks and N antiquarks to combine into M mesons, 

B baryons and B anti-baryons according to QCR is given by 

2N(N\f(M + 2B-l)\ oM _ 
XmB ~ (2N)\M\(B\y 3 6n ' M+3B (4) 

The average numbers of primary mesons M(N) and baryons B(N) are 

M(N) = E MX MB (N) 

M,B (5) 

B(N) = E BX MB (N) 

M,B 

Approximately, in the combination for iV > 3, M(N) and baryons B(N) can be well pa- 
rameterized as linear functions of quark number N, 

M(N) =aN + b 

(6) 

B(N) = G=2*JV- § 
where a = 0.66 and b = 0.56. But for N < 3, one has 

M(N) = N, B(N) = for N <3 (7) 

So that, the production ratio of baryon to meson is obtained from eqs. (jSJ) and 

Rb/M ~ 3(aN + b) (8) 
We see that SDQCM treats meson and baryon formation uniformly, and there is no extra 
ad hoc mechanism and free parameters for the baryon production. Here the B/M ratio is 
completely determined at a certain N, unlike in LSFM it is completely uncertain and has 
to be adjusted by free parameters, i.e, the ratio of diquark to quark qq/q. 



2. momentum distribution of primary hadrons in the qq system 



In order to give the momentum distribution of primary hadrons, each phenomenological 
model must have some inputs. For example, in LSFM, they use a symmetric longitudinal 
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fragmentation function 



/(z) oc (1 ^ a exp (-b^) (9) 



where a and b are two free parameters (and a is flavor dependent). In this paper, in order to 
give the momentum distribution of primary hadrons produced according to QPR and QCR, 
we simply adopt the widely used LPSA which is equivalent to the constant distribution of 
rapidity. Hence a primary hadron i is uniformly distributed in the rapidity axis, and then 
its rapidity can be written as 

Yi = z + £jr, o<&<i, (io) 

where & is a random number; Z and Y are two arguments, and can be determined by 
energy-momentum conservation in such a color singlet system 

J2E l = W, £P i4 = 0, (11) 

1=1 1=1 

where and Pn denote the energy and the longitudinal momentum of the zth primary 
hadron respectively, obtained by 

TP — m exp(Yj)+exp(-Yj) 

1 ~ Tl 2 (12) 



p _ ™ exp(Yi)-exp(-Yj) 
*IA — m Ti 



2 



where rriTi is given by 



m Ti = Y m\ + P Tl (13) 

— > 

where rrii is the mass of the zth primary hadron, and Pti obeys the distribution 

h 2 H 

/(Pti, ■ ■ ■ , Pth) oc J{exp{-^-)5(£PTi) (14) 

i=l ° i=l 

In this paper, we set cr = 0.2 GeV. Eq. (fl4"|) is just what LSFM uses. 

Note that LPSA or the constant rapidity distribution is rather naive, but it is convenient 
for us to study the correlations without introducing many complicating parameters. 



3. hadronization of a multi-parton state 



At the end of parton showering, a final multi-parton state will start to hadronize. To 
connect the final multi-parton state with SDQCM, we adopt a simple treatment assumed 
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in WCFM, i.e. before hadronization, each gluon at last splits into a q'q' pair, the q' and 
q' carry one half of the gluon momentum, and each of them forms a color singlet with 
their counterpart antiquark and quark in their neighbourhood, respectively. Now take the 
three parton state qqg as an example to illustrate the hadronization of a multi-parton state. 
Denote the 4-momenta for q, q, g as 

P l = (E q ,P q ), P 2 = (E g ,P g ), P 3 = (E q ,p q ). (15) 

Before hadronization, the gluon splits into a q'q' pair and the q' and q' carry one half of the 
gluon momentum, and the qqg system forms two color singlet subsystems qq' and q'q. The 
invariant masses of the subsystems are 



w m , = J & + f ) 2 = i(E q + fy - (p g + f ) 2 



(16) 



As was commonly argued by Sjostrand and Khoze recently [23], the confinement effects should 
lead to a subdivision of the full qq system into color singlet subsystems with screened inter- 
actions between these subsystems qq' and q'q. Hence SDQCM can be applied independently 
to each color singlet subsystem,^. e., we can apply the equations in the former two subsections 
to each subsystem, and obtain the momentum distribution for the primary hadrons in their 
own center-of-mass system. Then after Lorentz transformation, the momentum distribution 
of the primary hadrons in laboratory frame is given. This treatment can be extended to a 
general multi-parton state, for example, see Fig. HJ 

Obviously, when the emitting gluon is soft or collinear with the direction of q or q, qqg 
cannot be distinguished from qq and W qq > or W q > q is too small for hadronization. For the 
avoidance of these situations, a cut-off mass M min has to be introduced. Here M min is a 
free parameter in the perturbative phase. Its value and energy dependence is theoretically 
uncertain. The physical assumption is that M min is independent of energy [li [24 1. 



III. RESULTS AND DISCUSSION 

In this paper, we study the hadronization process at Z° factory by LSFM and SDQCM. 
We focus on investigating the properties of the heavy hadrons, e.g., the ratio of baryon to 
meson and the baryon-antibaryon correlations. For LSFM, we use the default parameter 
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O q • q 



FIG. 4: quark combination 



15| which fit 



values in PYTHIA, while for SDQCM, we adopt the parameters used in ref. 
data quite well. 

The final hadron multiplicity at high energy e + e~ reactions is always an important topic. 
We list the predictions of final hadron multiplicities at Z° factory = 91.2GeV) by 
LSFM and SDQCM in table I. The results of LSFM are obtained by running PYTHIA6.4 



with the default parameters. The experimental data are from ref. [25]. It is found that 
the predictions of LSFM and those of SDQCM are consistent with most of the experimental 
data. However, the multiplicity of the heavy baryons, eg., H&, £&, are still not measured 
at Z° energy, and the corresponding theoretical predictions obtained by LSFM and SDQCM 
are quite different. Therefore it is important to measure their production rates at Z° pole in 
order to discriminate different hadronization mechanisms. At the Z° energy, we study the 
relation between the integrated luminosity and the producing number of the heavy baryons. 
The corresponding results are displayed in Figj5l It is obvious that for the LEP I, the 
production number for E b predicted by LSFM(SDQCM) is about 10 3 , and that for Q b is 
several tens (hundreds). In order to study the heavy baryon production mechanism with 
enough precision, the integrated luminosity should be increased as large as possible. For 
example, if the integrated luminosity of Z° factory can reach about 10 4 p6 _1 , the produced 
number of Qb may reach several thousands(ten thousands) according to the prediction of 
LSFM(SDQCM), which makes it possible to study the properties of the heavy baryons and 
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to test the hadronization models with higher precision. 
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FIG. 5: The relation between integrated luminosity and the production number of heavy bayons 
at Z° energy. The lower solid line stands for the total luminosity of LEP I from 1993 to 1995 at 
energies around the Z° 261 ] . 



We also show the results for the momentum distribution of 7r ± , K ± and pp in Figj6j 
The results of LSFM are more consistent than SDQCM. This is easy to understand, since 
SDQCM does not include any complicated inputs for the hadron momentum distributions 
to be improved in the future. 

The following focuses on studying the baryon to meson ratio and BB correlations which 
reflect the hadronization mechanism more directly. In LSFM, the baryon to meson ratios 
can be tuned by the free parameters, e.g., qq/q, BMB j[BB + BMB). SDQCM describes 
the baryon and meson production in the uniform scheme, so that these ratios are obtained 
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FIG. 6: The momenta spectrum of ir^, and pp at Z° factory. Solid(dashed) line is for 
LSFM(SDQCM), and here Pbeam = VS/2. The experimental data is from ref. 27^ . 

naturally. Some of these ratios are listed in table II at Z° factory. One can notice that both 
LSFM and SDQCM can explain the present data. In e + e~ annihilations, most baryons are 
produced directly, even if they are the decay products. They carry the main properties of 
their mother particles. 

Studying the properties of baryons, especially the BB flavor correlations, is helpful to 
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reveal the hadronic mechanism. Here the BB flavor correlation strength is defined as 

R bS = N - Upair 



n B + n§ 



where n pa i r is the number of the BB pairs and Nb(N§) is the baryon(antibaryon) number. 
N = 2(1) in the case of B = B(B ^ B) (27)]. The predictions of baryon antibaryon flavor 
correlation strength by SDQCM and LSFM are listed in table III together with the corre- 
sponding OPAL data[27|. The understanding of hadronization mechanism related to the 
heavy hadron production requires more precise measurements at the future Z° factory. 



IV. SUMMARY AND OUTLOOK 

In the e + e~ — > 7*/^*° — > h! process, the final particles have no relations with the structure 
of the initial particles, so it is convenient to study the hadronization mechanism. Above all, 
such studies help to investigate the multi-production processes in other high energy reac- 
tions^.^., at RHIC energy). In this paper, we report the predictions especially related to 
heavy hadrons, for example, baryon to meson ratio and BB flavor correlations obtained by 
LSFM and SDQCM in the hadronic Z° decay. These results show that the future exper- 
iments at Z° factory to test hadronization mechanism, especially rare hadronized events, 
such as the doubly heavy baryon production, are significant. 

Though many studies have been made for the e + e~ — > h's process, still a lot of unclear 
but important topics are worth studying, such as the color connection effects 2s|. The color 
connection among final partons plays a crucial role in the surface between PQCD phase 
and the hadronization one (see Fig. [TJ. Different choices of the color connections of the 
final multi-parton system leads to different hadronization results. Recently, the doubly 
heavy baryon(e.g., S cc ) production has been investigated in different reactions [29]. In e + e~ 
annihilation, the doubly heavy baryon production reveals a special kind of color connection 
of the final parton system q\q\qi<l2 + ng(n < 0). This kind of color connection has not been 



considered in the popular hadronization models. In ref. |30j], the hadronization effects of this 
kind of color connection are investigated in an extremely limited case of e + e~ — > ccqq — > h's 
process at B factory energies. There cq and cq form diquark-antidiquark pair, which is 
treated as a color singlet string and fragments into hadrons in the framework of LSFM. 
However due to the phase space limitation, the hadronization effect induced by diquark pair 
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fragmentation is small. Fortunately, in the 
as e + e~ — > Z° — > E cc ccX will be produced 



uture Z° factory, a large number of events such 



31|, so that the correlated physical observables 



can be measured with higher statistics. As a result, the hadronization mechanism and the 
related non-perturbative nature can be well studied at Z° pole energy. 
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TABLE I: Results for average hadron multiplicities at Z factory. The experimental data are from 
ref. [25|]. 



Particle 


EXP DATA 


LSFM 


SDQCM 


7T+ 


17.02 ±0.19 


17.125 


17.766 


7T° 


9.42 ±0.32 


9.696 


9.633 


K+ 


2.228 ±0.059 


2.312 


2.145 


K° 


2.049 ±0.026 


2.079 


1.758 


7] 


1.049 ±0.080 


1.013 


0.787 


D+ 


0.175 ±0.016 


0.166 


0.217 


D° 


0.454 ±0.030 


0.495 


0.476 


B+ 


0.178 ±0.006 


0.174 


0.182 


B° 


0.165 ±0.026 


0.173 


0.182 


B° s 


0.057 ±0.013 


0.052 


0.053 




1.016 ±0.026 


1.369 


1.809 


P° 


1.231 ±0.065 


1.524 


1.842 


K*+ 


0.715 ±0.059 


1.112 


1.043 


K*° 


0.738 ±0.024 


1.107 


1.012 


D*+ 


0.1937 ±0.0057 


0.2399 


0.2547 


D*° 




0.2407 


0.2603 


P 


1.050 ±0.032 


1.223 


0.928 


A 


0.3915 ±0.0065 


0.3930 


0.4368 


E° 


0.076 ±0.011 


0.075 


0.135 


E - 


0.081 ±0.010 


0.069 


0.114 


S+ 


0.107 ±0.011 


0.074 


0.122 


A c + 


0.078 ±0.017 


0.060 


0.076 


Ag 


0.031 ±0.016 


0.035 


0.044 


E° 




0.0017 


0.0073 






0.0019 


0.0102 






0.0024 


0.0065 






0.00006 


0.0008 
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TABLE II: Results for baryon to meson ratios at Z factory. The experimental data are from ref. 
25]. 





EXP DATA 


LSFM 


SDQCM 


A? 

w 


0.188 ±0.101 


0.201 


0.239 


B+ 


0.174 ±0.090 


0.200 


0.240 


At 
~W 


0.172 ±0.039 


0.121 


0.160 


At 
IF 


0.446 ±0.105 


0.360 


0.351 



TABLE III: Predictions of baryon antibaryon flavor correlations by SDQCM and LSFM, and the 
experimental results are from OPAL collaboration|27l] . 

EXP DATA LSFM SDQCM 

AA 0.49 ±0.06 0.38 0.48 

E~E+ 0.04 ±0.06 0.14 0.15 

Z-A + Z+A 0.463 ±0.099 0.510 0.538 

A° b A° b 0.08 0.12 
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